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ABSTRACT 

The spin temperature of neutral hydrogen, which determines the 21 cm optical depth and 
brightness temperature, is set by the competition between radiative and collisional processes. 
In the high-redshift intergalactic medium, the dominant collisions are typically those between 
hydrogen atoms. However, collisions with electrons couple much more efficiently to the spin 
state of hydrogen than do collisions with other hydrogen atoms and thus become impor- 
tant once the ionized fraction exceeds 1%. Here we compute the rate at which electron- 
hydrogen collisions change the hydrogen spin. Previous calculations included only S-wave 
scattering and ignored resonances near the n — 2 threshold. We provide accurate results, in- 
cluding all partial wave terms through the F-wave, for the de-excitation rate at temperatures 
Tk < 1 .5 X 10"* K; beyond that point, excitation to n > 2 hydrogen levels becomes significant. 
Accurate electron-hydrogen collision rates at higher temperatures are not necessary, because 
collisional excitation in this regime inevitably produces Lya photons, which in turn dominate 
spin exchange when Tk ^ 6200 K even in the absence of radiative sources. Our rates differ 
from previous calculations by several percent over the temperature range of interest. We also 
consider some simple astrophysical examples where our spin de-excitation rates are useful. 

Key words: atomic processes - scattering - diffuse radiation 



1 INTRODUCTION 

The 21 cm transition is potentially a powerful probe of the pre- 
reionization intergalactic medium (IGM) because of the eno rmous 
amount of neutral hydrogen in the Universe at that time ('Fielcf 
11958; Scott & Rees 1990; Madau et al. 1997). It can teach us about 
reionization, the formation of the first structures and the first 

I alaxies, and even the "dark ages" before these objects formed 
Furlanetto et all2006l and references therein). It is therefore cru- 
cial to understand the fundamental physics underlying the 21 cm 
transition. One critical aspect is the spin temperature, which is de- 
termined by the competition between the scattering of cosmic mi- 
crowave back ground (C MB) photons, the scattering of Lya pho- 
tons I Wouthuvsei] l 1 9521; fpield 1958), and collisions. When CMB 
scattering dominates, the IGM remains invisible because the spin 
temperature approaches that of the CMB (which is used as a back- 
light). 

Before star formation commences, collisions are the only 
way to break this degeneracy. The total coupling rate is deter- 
mined by collisions both with other hydrogen atoms and with elec- 
trons. At the low resid ual electron fractio n expected after cosmo- 
logical recombination ^SeageretalJll999^ ■ H-H collisions domi- 
nate. Spin exchange in such interactions has received a great deal 
of attention over the years (.Purcell & Field. ■195& ; .Smith .1966 ; 
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I Allison & Dalgamol 1 19691; IZygelmaiJ I200I iHirata & Sieurdsonl 
2006J). Such collisions suffice to couple the spin temperature Ts 
to the kinetic temperature Tk at high densities, when z^50 
{ Loeb & Zaldar riaga 2004). However, the spin exchange rate in 
these collisions falls dramatically at Tk ^ 100 K; coupled with 
the declining density, such interactions become ineffective at 
lower redshifts. H-e^ collisions have received much less at- 
tention. However, at fixed temperature, free electrons move at 
much higher velocities than do hydrogen atoms, and are thus 
much more efficient at spin exchange on a per-particle basis. If 
the mean ionized fraction is Xi > 0.01, they can play a signifi- 
cant role in the coupling. This situation can occur if, for exam- 
ple. X-rays from the first quasars or sta r-forming galaxies cre- 
ate a warm, partially-ionized medium tOhll200lHVenkatesan et alJ 
[ 200 it iRicotti & Ostrikedl2004 iRicotti et alJl2005t lNussedl200a 
iKuhlen & Madau"200^ 'Kuhlen et al.'200d^, or if the decay or an- 
nihilation of exotic particles deposits a significant a mount of en- 
ergy into the IGM (Chen & Kamionkowski 2004; Pi erDaolill2004 
iShchekinov & Vasiliev«200& ; .Furlanetto et alj.2006) . 

The spin de-excitation rate in H-e~ collisions, which we will 
denote kiq, obviously depends on the detailed scattering proper- 
ties of these collisions. Over the past several decades, this inter- 
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trons have energies below the threshold to excite the n = 2 level, 
and the dominant mechanism of spin de-excitation is electron- 
electron spin exchange. Although this regime has been particu- 
larly well-studied, it has only been solved in full nonrelativi stic de- 
tail re latively recently. The latest spin exchange calculation JSmithI 
Il96^ use d the scattering phase shifts computed numerically by 
ISchw artz ( 1961). While accurate, these included only S-wave scat- 
tering and neglected both higher-order partial waves and H~ reso- 
nances near the n = 2 threshold (which were not resolved by the 
numerical algorithm). Both of these issues have been solve d over 
the intervening four decades (see, e.g., Wang & Callawavl ll993l 
|l994 for a recent calculation). Our goal is to recalculate k\o in 
light of these newer phase shifts. 

Higher-order partial waves become significant at 
Tk ^ 1000 K; if X-ray heating in the early Universe is strong, 
then these temperatures can be achieved easily (Furlanetto 2006). 
We therefore examine the high-temperature limit in some detail. 
Our calculation breaks down once coUisional electronic excitation 
(and ionization) become important. We will show that, above 
Tk ~ 1.5 X 10* K, direct spin de-excitation by scattering above 
the n = 2 threshold cannot be ignored. Although such interactions 
are difficult to model - because of the wide array of possible 
transitions - in practice this regime is relatively unimportant. In the 
low-density IGM, coUisional excitation to the 2p and higher levels 
is followed rapidly by radiative de-excitation, producing a Lya 
background. We show that this background inevitably dominates 
the spin temperature coupling at Tk > 6200 K (assuming a 
thermal distribution of electrons). The crossover can occur at 
smaller temperatures if a non-thermal po pulation of fast electrons 
exists (produced, for example, by X-ravs: IChen & Miralda-Escudd 
12006; Chuzhov et a], 2006). 

The remainder of this paper is organized as follows. In <J2| 
we briefly review the 21 cm transition. Our main results are con- 
tained in S|3| where we examine various mechanisms for spin de- 
excitation, consider the high-temperature limit, and calculate the 
spin de-excitation rates for H-e^ collisions. We present some sim- 
ple astrophysical applications in Sj4| and we conclude in J5| 

In our numerical calculations, we assume a cosmology with 

= 0.26, 9.K = 0.74, ilt = 0.044, H = WOh km s"^ Mpc"^ 
(with h = 0.74), n — 0. 95, and as = 0-8, consistent with the 
most recent measurements jSpereel et all2006h . although we have 
increased ag from the best-fit WMAP value in order to improve 
agreement with weak-lensing data. 



The spin temperature Ts is determined by competition be- 
tween scattering of CMB photons, scattering of UV photons 
( Wou thuysen 1952; Field .1958i) . and collisions >Purcell & Fiei3 
Il956l) . In equilibrium,^ 



rp—l T^ + XqTj^ + XaTf. 

a. 



(2) 



Here Xc is the total coUisional coupling coefficient, including both 
H-H and H~e~ collisions. We will break the total coupling into two 
components, x™ and x'^, with obvious meanings. The coefficient 
from H-e^ collisions is 



cH UeKloTi, 



(3) 



where rie is the local electron density, T, = hv2i/kB ~ 0.068 K, 
V21 is the frequency of the 21 cm line, and Aio is the Einstein- 
A coefficient for that transition. The last part of equation J2} de- 
scribes the Wouthuysen-Field effect, in which the absorption and 
re-emission of Ly a photons mixes the hyperfine states. The cou- 
pling coefficient is iChen & Mi^alda-Escud^2004^ 



= 1.81 X 10"(l + z)"^5'aJc, 



(4) 



where Sa < 1 describes the detailed atomic physics of the scat- 
tering process and Ja is the background flux at the Lya fre- 



quency (ignoring scattering) in units of cm s 



-iHz- 



the 



Wouthuysen-Field effect becomes efficient when there is ~ 0.1 
photon per baryon near this frequency. It couples Ts to an ef- 
fective color temperature Tc (in most circumstances, Tc ~ Tk', 
[Field 1959) . Several estimates of Sg and Tc exist in the literature 
Ichen & Miralda-Escude 2004; Hiratj |2006l: IChuzhov & Shapiro! 
[2005^ .Furlanetto & PritcharA2006) . 



3 ELECTRON-HYDROGEN COLLISIONS 

3.1 Mechanisms of coUisional spin de-excitation 

In principle there are five distinct coUisional mechanisms that can 
cause spin de-excitation. The existing literature has addressed only 
one of these, electron spin exchange. As we intend to increase the 
accuracy to which kio is known, we will first reconsider the relative 
magnitudes of each mechanism. We will then calculate spin de- 
excitation cross sections and compare them to previous results. 



2 THE 21 CM TRANSITION 

We review the relevant characte ristics of the 21 c m transition here; 
we refer the interested reader to lFurlanetto et alj i2006 ) for a more 
comprehensive discussion. The 21 cm brightness temperature (rel- 
ative to the CMB) of a patch of the IGM is 



Sn = 27xHi(l + cS) 

'Ts-1 
Ts 



nth^ \ / 0.15 l + 2\i/2 
0.023 ) [a^ ~iO 

H{z)/{l + z) 



mK, 



(1) 



where 5 is the fractional overdensity, shi is the neutral fraction, 
Xi = 1 — xm is the ionized fraction, Ts is the spin temperature, T-y 
is the CMB temperature, and dwy /dry is the gradient of the proper 
velocity along the line of sight. The last factor accounts for redshift- 
space distortions iBharadwai & Ali.2004 : .Barkana & Loeb.2005.) . 



3.1.1 Electron-electron spin exchange 

Electron-electron spin exchange refers to the process by which the 
spin of the incoming electron is exchanged with that of the atomic 
electron. One such collision could be represented schematically by 



q(s) + P{a)cx{p) -> /3(s) + a(a)a(p), 



(5) 



where a and (5 refer to the Pauli spin states of spin-i particles, s 
and a refer to the scattering and atomic electrons, respectively, and 
p refers to the proton. In this collision, a singlet hydrogen atom is 



^ This is an excellent approximation throughout cosmic history. The worst 
case occurs when interactions with CMB photons dominate; the relevant 
timescale is then ^ (Bio/^)^^, where _Bio is the Einstein absorption 
coefficient and is the CMB intensity at the 2 1 cm transition. This yields 
i'y/iH ~ 4 X 10~''(1 + z)^/^, where is the Hubble time. If colHsions 
or Lyo absorption dominate, the timescale will obviously be even smaller. 
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converted into a triplet hydrogen atom by spin exchange. Includ- 
ing all possible spin permutations of this collision leads to an ex- 
pression for the spin de-excitation cros s section. (For a more de- 
tailed discussion of spin exchange see iMott & Mas'sevlll965l and 
ICondon & ShorflevI 19631 ) 

Ignoring all relativistic interactions and non-central forces, 
the spin-exchange-mediated spin de-excitation cross section can be 
calculated relatively simply. (Interactions which violate these as- 
sumptions will be discussed below in i|3.1.2l and i|3.1.3l ) We be- 
gin with the time-independent Schrodinger equation, which may 
be written 



//*(x,y) = i5*(x,y) 
where, in atomic units. 



^2 2 2 



(6) 



(7) 



m \y\ |x-y| 

is the total wavefunction, and x and y label the positions of the 
two electrons. The Pauli exclusion principle demands that 

*(x,y) = ±*(y,x), (8) 

where the spatially symmetric and antisymmetric versions of the 
equation respectively refer to singlet and triplet scattering. We treat 
the scattering as electronically elastic (see il3.1.5t . so the solution 
has the asymptotic form of 

*~e*-^ + /k(e)^, (9) 



where z is the direction of the incoming electron, k is its momen- 
tum, r is the magnitude of the vector r separating the two electrons, 
and fk{d) is the scattering amplitude at angle 6. Thus, the differ- 
ential scattering cross section into solid angle dQ is da/dQ — 

We solve for the elastic scattering amplitude by expanding the 
incoming electron into partial waves of different orbital angular 
momentum (e.g., Mott & Massev 1965), yielding 



As a three-body interaction, no analytic solution exists, although 
one can be found in the restricted problem of zero total and or- 
bital angular momentum (Temkin 1962; Poet 1978). It has been 
solved through variational principles ( Schwartz 1961; Shimamup] 
I197J; iDas & Rudgelll976l ; iReeister & Po&.I975; .Callawav..l97^ 
the close coupling formalism (in which the total wavefunction 
is exp anded in a basis constructed from hydrogen eigenstates; 
iBurke & Schey. ,,1962a.), the related convergent close coupling 
formalism iBravetalJl2002h. a nd through fully numeric tech - 
niques jWang & Callawavlll993l Il994 IShertzer & Boterdll994) . 
All methods now agree on to several significant figures in 
the range 0.1 < A: < 0.8; furthermore, these results agree with 
laboratory experi ments on spin-excha nge effects in higher-energy 
H-e~ collisions jFletcher et alJll985h . We will use the results of 
IWang & CaIlawa^^ j 1994^ 7 who present the L = 0-3 singlet and 
triplet phase shifts. Neither the D- nor the F-wave term contributes 
significantly to a^c, so truncating the sum at i = 3 is sufficient for 
our purposes. 

The cross section at zero energy requires more subtle anal- 
ysis, especially because the behavior at fc < 0.1 (or i5 < 0.01 
Ry) is crucial for Tk ^ 1000 K. In the low-energy limit, centrifu- 
gal barriers make the L > terms vanish, but the 5-wave scattering 
term remains finite. The solution is typically presented in terms of 
the scattering lengths as.t, defined so that tan^^'* ~ kas,t- Then, 
again in units of Aq, jSeatonll957h 



lim Use = -r(at 
k^o 4 



(14) 



Scattering lengths are difficult to compute because the effective po- 
tential seen by the electron at zero energy dies off rather s lowly. 
The most recent calculations appear to be bv lSchwartzlil96ll) . who 
found that = 5.965 ± 0.003 and at = 1.7686 ± 0.0002. 

We present detailed numerical results for the spin-exchange- 
mediated spin de-excitation cross section in il3.3l below. For com- 
parison with other mechanisms, the magnitude of these cross sec- 
tions in atomic units is of order unity. 



(10) 



where Pl is a Legendre polynomial of order L, S'^ is the phase 
shift for angular momentum L in the singlet (s) and triplet (t) spin 
states, and k is in units of the Bohr radius ao (in this system ener- 
gies are naturally expressed in Rydbergs). To obtain the total cross 
section, we simply average over the initial spin states of the hydro- 
gen atom, so (in units of ag) 



oc 

'Ttot(fc) = ^ ^(2L + 1) (i sin' Sl + ^ sin' Sl 



(11) 



The spin exchange cross section a^c is more complicated, be- 
cause it is a coherent s um of the scattering amplitudes (.Fieldkl958i ; 
iBurke & Schevll962bt) : 



so that (again in units of gq) 

00 

aso{k) = ^ ^(2L + 1) sin'(5i - SI 



(12) 



(13) 



Thus, calculating the spin exchange cross section only re- 
quires knowledge of the phase shifts for elastic scattering. This 
is a straightforward, though computationally challenging, problem. 



3.1.2 Electron-electron interactions 

In addition to spin exchange, the incoming electron can interact 
with the inherent magnetic dipole moment of the atomic elec- 
tron. (Note that spin-i systems have only monopole and dipole 
moments, so higher order multipoles need not be considered; 
iBethe & BacheJl93d ) 

The incoming electron interacts with the magnetic moment of 
the atomic electron by inducing a torque through a magnetic field. 
The interaction Hamiltonian, i/int-c, is given by 



Hi 



lit — c — 



Me 



Me 



(15) 



where L is the orbital angular momentum of the incoming electron, 
/x^ is the magnetic moment of the electron, given by 

gech 



Me = 



(16) 



4mc ' 

is the Lande factor for the electron, the other symbols have their 
usual meanings, and we have ignored any effect of the atomic nu- 
cleus I Jackson 1999). The first three terms arise from the intrinsic 
magnetic moment of the incoming electron, while the final term is 
generated by the motion of the incoming electron and its associated 
charge. A full solution of this interaction would require relativistic 
quantum scattering theory with non-central forces and is beyond 
the scope of this study. 
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Fortunately, relativistic effects such as spin-spi n and 
spin-orbit coupling scale as jc? iSerestetsk ii et alj Il982t 
ICjordrnan & Krivchenkov 19931) . As an upper limit, an electron 
with sufficient energy to ionize the hydrogen atom (1 Ry) has ve- 



locity v = 2 X 10 cm s 



5 X 10"^ (We will 



show below that this velocity corresponds to a temperature above 
those at which our model is valid, so it suffices for a worst case 
estimate.) Compared with the cross section of order unity due to 
the spin-exchange mechanism, this is ne gligible. The good ag ree- 
ment between the experimental results of lFletcher et al and 
the cross sections calculated from spin-exchange alone further sup- 
ports our subsequent neglect of these interactions. 



3.13 Electron-proton interactions 

The incoming electron can also interact with the spin of the atomic 
proton. Again, this can occur in two different ways - through the 
spin of the electron and through the magnetic field generated by its 
motion. The interaction Hamiltonian, Mnt-p, is given by 



87r(/i^_/ip)(5(r) 
3 

(Me ■ Mp) 



+ 



,(r-Me)(r-Mp) 



where {i^ is the magnetic moment of the proton, given by 



(IV) 



(18) 



rrip is the mass of the proton, is the Lande factor for the proton, 
and the other symbols are defined as in equation <15t . 

The same arguments used above for the electron-electron in- 
teraction apply to the electron-proton interaction. In this case, how- 
ever, the magnetic moment of the nucleus is smaller than that of the 
electron by the ratio of their masses, ~ 1/1836, so the cross section 
is further suppressed to order 10^*, and neglecting this interaction 
will not alter our results. 



3.1.4 Transient complex formation 

In the region k > 0.8, a series of resonances corresponding to 
quasistable doubly excited states of occur. These resonances 
have two effects on our calculations. First, the elastic scattering 
phase shift undergoes an abrupt change by vr at a resonance. The 
resonant phase shift can be fit to the usual Breit-Wigner form 



tan((5 — 



r/2 

Eb. — E ' 



(19) 



where St is the smoothly varying background phase shift, Er is 
the center of the resonance, and F is its width. In order to estimate 
the effects of this structure on the spin-exchange-mediated spin de- 
excitation cross section, we include the lowest order ^S, '^P, a nd 
D resonances according to the fits of lWang & CalTawavlJl994 . 
Second, as the system passes through these resonances, it tran- 
siently becomes H~ . For a full treatment, we would need to con- 
sider the evolution of the spin states during the time spent as H~ 
and the subsequent autodetachment reaction. Two facts lead us to 
ignore this possibility. First, the resonances are broad. The narrow- 
est is the state, with a width of 2 x 10"* Ry JWang & Callawavl 
Il994) : the corresponding natural lifetime of 12 fs is much shorter 
than typical nuclear spin relaxation times. Second, the resonances 
cover less than half of a percent of the energy range considered. To 




0.4 

E (Rydbergs) 

Figure 1. Spin-exchange cross section (in units of ttOq) for H-e~ colli- 
sions, as a function of the energy of the incident electron. The solid curve 
shows (Tse including partial waves with L < 3, while the dashed curve 
includes only the L = term. We include the three lowest energy H~ 
resonances. 



the level of accuracy we seek, even complete spin de-excitation at 
these resonances will not change our results. 



3.1.5 Electronic excitation of hydrogen 

At the low temperatures relevant to the high-redshift IGM, nearly 
all collisions occur below the n = 2 excitation threshold. However, 
once this threshold is reached, excitation to higher-n levels quickly 
becomes energetically feasible. Past this point, the plethora of pos- 
sible transitions makes it difficult to compute the H-e" collisional 
spin-exchange rates from first principles. We can, however, esti- 
mate the temperatures at which such excitations can affect the de- 
excitation rate significantly. Detailed numerical results are given 
in i|3.3l the conclusion is that a one percent effect is reached at 
Tk ~ 1.5 X 10* K. 

In practice, kw quickly becomes irrelevant once the n = 
2 excitation threshold is reached, because the Lya background 
generated by collisional excitations and subsequent radiative de- 
excitations will then dominate the spin de-excitation process. The 
production rate of Lya photons (in units of photons per volume per 
second) is ~ nuineKa, where Ka is the rate coefficient for exci- 
tations that eventually produce Lya photons and nm is the density 
of hydrogen atoms. For a simple estimate, we set Kc, — K2p, the 
rate of direct excitations to the 2p level. This is a lower limit be- 
cause of cascades from higher levels: note that although roughly 
one-third of excitations to higher p states produce Lya photons 
(Hirata 2006; Pritchard & Furlanetto 2006b a), the cross sections 
for these excitations are significantly smaller than the cross sec- 
tion for the ls-2p transition, so the approximation is justified. (For 
example, is nearly fifty times smaller than K2p-) Under this ap- 
proximation, the background flux at the Lya resonance (assuming 
neutral helium) is 
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Figure 2. {a): Rate coefficients. The tiiicJc solid line shows kio, the spin de-excitation rate from elastic collisions below the n = 2 thi'eshold. The thin solid 
line shows the same quantity if we only include the L = partial wave. The long-dashed and dot-dashed curves show k^^2 and k „>^, respectively. The 
short-dashed curve shows K2p- For comparison, the dotted curve shows ■ (b): The solid curve shows the ratio of the lLiszll l200ll) fit for kiq to the exact 
result. The dashed curve shows the ratio of the L = solution to the full result. 



(z)K2p 



(20) 



where Va is the frequency of the Lya transition. In our cosmology, 
the coupling coefficient is then 

+ 



/ ^2p 

V2.5 X 10-15 cm^ s^ 
a:Hi(l - xm 
0.25 



(21) 



Using the cross sections calculated below, we find that K2p is suffi- 
ciently large for Lya scattering to dominate at Tk ~ 6200 K; see 
i|3.3l for a more detailed discussion. 

Overall, then, the effect of electronic excitations of hydro- 
gen is to limit the applicability of the present calculations at high 
temperatures (above Tk ~ 1-5 x 10* K). Between that limit and 
Tk ~ 6200 K, the effects of Lya photons must also be included. 

3.2 The cross section 

Having considered all the possible mechanisms for spin de- 
excitation, we have found that the spin-exchange mechanism does 
indeed dominate under astrophysical conditions, at least to the frac- 
tion of a percent level. There are possibly some effects of H" res- 
onances in the region around fc = 0.8 Ry, but otherwise our results 
are valid to better than the percent level up to Tk ~ 1-5 x 10* K. 

In this section we proceed to calculate actual spin de- 
excitation cross sections from the spin exchange rates. Our inclu- 
sion of higher-order partial waves is the principal improvement 
over previous calculations. Field ( 1958) used the approximate S- 
wave phase shifts of Massev & Moiseiwitch (1951), which differ 
from the true values by several percent over this energy range. Most 
significantl y, their scattering lengths differ by ~ 25% from the cor- 
rect values. Smith ( 1966) used the (near-exact) S-wave phase shifts 
of lSch wartz ( 1961) but did not include the L > Q terms. 



FigureQshows cjac as a function of collision energy. The solid 
curve includes the L — 0-3 terms, while the dashed curve includes 
only L = a nd thus shows nearly the same cross section used by 
ISmithl il96d) . Clearly the higher-order partial waves become im- 
portant at fc>0.1, where non-zero angular momentum collisions 
become common. In fact, the P-wave makes up the vast majority 
of the difference (had we included only L — and 1, the solid 
curve would appear the same except for the highest-energy reso- 
nance). Note the resonance structure at fc > 0.8; fortunately, al- 
though Use changes rapidly in this regime, the resonances constitute 
only a small part (F < 3 x 10^'' Ry) of the energy range. 



3.3 The rate coefficient 

The H-e~ coUisional spin de-excitation rate coefficient is 



ISksTK _ 



(22) 



where the prefactor is the mean collision velocity, M ~ me is 
the reduced mass of the II~e~ system, and the thermally-averaged 
cross section is 



1 



(fcsTK)2 



dE asciE)Ee' 



-E/kgTK 



(23) 



Similarly, we can define k„<2, the rate coefficient for all collisions 
below the n = 2 threshold, with the replacement asa — > o"tot- 

The heavy solid line in Figure|2jj shows kio if we include only 
collisions below the n — 2 threshold. (In other words, we truncate 
the integral in eQ. I23l at 10.2 eV; this restriction causes the decline 
at Tk ^ 40, 000 K.) At low temperatures, it is nearly proportional 

1/2 

to T^' (which comes from the velocity factor): scattering is dom- 
inated by electrons near zero energy, where the amplitude is es- 
sentially set by the scattering lengths a^^t. kiq flattens out at higher 
temperatures because aae declines with energy (see Fig.0. We note 
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Table 1. Electron-hydrogen spin de-excitation rates 





KlO (10"^ cm^ s~i) 


Tk (K) 


Kio (10~^cm3 s~i) 


1 
1 




lUUU 


j.yz 


2 


0.337 


2000 


7.15 


5 


0.530 


3000 


7.71 


10 


0.746 


5000 


8.17 


20 


1.05 


7000 


8.32 


50 


1.63 


10,000 


8.37 


100 


2.26 


15,000 


8.29 


200 


3.11 


20,000 


8.11 


500 


4.59 







that the resonances in FigureQare so narrow that they never affect 
Kio by more than 1%. In contrast, the dashed curve shows the to- 
tal rate k„<2 for all collisions below the n = 2 threshold. Only 
~ 10% of collisions actually result in a net change of hydrogen 
spin. 

For comparison, the dotted curve in Figure|2jz shows the corre- 
sponding de-excitation rate co efficient for H-H collisions. The data 
at Tk ^ 300 K are taken from Zvselman ' 12005'). At higher temper- 
atures, we use cross sections from K. Sigurdson (private commu- 
nication, also tabulated in Furlanetto et al. 2006). We do not show 
at Tk > 10^ K because Sigurdson did not include excita- 
tion to n > 2 states from H-H collisions. At high temperatures, 
kio/k™ ~ mp/me (comparable to the relative velocities of 
the scattering species). The hydrogen cross section is much smaller 
at Tk ^100 K because of an accidental cancellation in the S-w ave 
cross sections iZvgelmarj2'00llSigurdson & Furlanettol200dl . 

The thin solid curve shows the cross section if we inclu de only 
the L = term (this is essentially identical to the result of lSmithI 
Il966ft : we also show the ratio between this version and our result 
by the dashed curve in Figure|2j'- As expected, higher-order partial 
waves only affect kio at Tk ^ 1000 K, where a substantial fraction 
of the electrons have large enough momenta for the L — 1 term to 
become significant. 

Calculation s in the literature often use the fit to kio proposed 
byjLiszj |200lt, which is based on the high-temperature results 
of ISmithI (1966). The solid line in Figure |2j' compares this func- 
tion to the exact value. At low temperatures, it exceeds our re- 
sults by several percent, and it systematically underestimates kio at 
T > 3000 K. We therefore recommend interpolating the exact rate 
coefficients; to this end. Table Q presents our results in numerical 
form. 

Unfortunately, there are no accurate calculations of a^a above 
the n = 2 threshold; typically only the total cross sections for 
elastic scattering, excitation, and ionization are presented (from 
which the spin-exchange cross section - a coherent sum of the 
singlet and triplet scattering amplitudes - cannot readily be ex- 
tracted). However, we can at least estimate the total rate (includ- 
ing elastic scattering, excitation, and ionization) for interactions 
above threshold (which we will call k„>2); comparing that to k„<2 
provides an estimate of the temperature at which our calculation 
breaks d own. To compute o ' „>2, w e interpolate the cross sections 
given bv lWang & Callawavl Il994l) for energies between the n = 2 
and n = 3 thresholds (which include elastic collisions as well 
as excitations to the 2s and 2p levels) and the total cross sec- 
tions for E > 13.6 eV from the Convergent Close Coupling on- 



10-11 , , , 




4x103 6x103 8x10= 10" 



Tk (K) 

Figure 3. Rate coefficient for H-e^ collisional excitation to the 2p level. 



line database'^ (see lBrav et all20o3 and references therein), which 
includes elastic scattering, excitation to all levels through 4/, and 
ionization. 

The dot-dashed curve in the left panel of Figure |2| shows the 
resulting rate coefficient k,„>2- It increases rapidly at Tk > 10* K, 
when the tail of the electron velocity distribution begins to populate 
theiJ> 10.2 eV space. We find fi:„>2/K„<2 ~ 0.01, 0.1, and 1 at 
Tk « (1.5, 2.3, 6.0) x lO" K. As mentioned above, we therefore 
expect that our mo is accurate to better than one percent at Tk < 
1.5 X 10* K. 

The short-dashed curve in Figure|2lj shows K2p (computed in 
the same way as the total cross section), which we used in ^3.1.51 
above to calculate the temperature at which Lya scattering be- 
comes important; K2p is typically several percent of /t„>2. Obvi- 
ously it is an extremely steep function of temperature, so we show 
a closeup of the regime of interest in Figure |3| We find that K2p 
reaches the level required by equation <21> at Tk ~ 6200 K. At 
Tk ^ 5000 K, the Lya background can be neglected in most appli- 
cations; on the other hand, by Tk ~ 8000 K, it will entirely dom- 
inate. In detail, the total Lya production rate will be slightly larger 
than s hown here because of radiative cascades from higher levels 
i Pritchard & Furlanetto 2006a,) . However, the extremely steep de- 
pendence of this rate coefficient on temperature suggests that in 
practice such minor corrections will not be very important. We also 
note again that Figure |3| assumes a thermal distribution of elec- 
trons; if a nonthermal population of fast electrons exists (as is the 
case if X-rays perm eate the IGM; Ichen & Miralda-Escuddl2006l 
IChuzhov et alJ2006h . they can collisionally excite higher levels and 
so produce Lya photons even if the mean IGM temperature is much 
smaller. Thus Lya production is probably always somewhat im- 
portant, but its details depend on the extremely uncertain radiation 
backgrounds (■Sethi,.2005. ; .Fiirlanettaii2006.) . We will therefore not 
address this possibility here. 



2 See |http://atoin.murdoch.edu.au/CCC-WWW/iiidex.html| 
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Figure 4. 21 cm brightness temperature fluctuations at 2 = 12 for Tk = 
10^ K. The thin curves take Xi = 0, 0.01, 0.05, 0.1, and 0.2, from bot- 
tom to top, with XoL = 0. The upper thick curve shows the fluctuation 
amplitude if Tg = lO'^ K, while the lower thick curve shows the fluctua- 
tion amplitude if the temperature and ionization fraction are those expected 
in an adiabatically cooling Universe. 



4 ASTROPHYSICAL APPLICATIONS 

Although Kio is much larger than the corresponding quantity for 
H-H collisions, H-e^ collisions are often ignored in calculating 
the spin temperature of the 21 cm line in the high-redshift IGM. At 
the small residual ionized fractions (xi ~ 2 x 1 0^"^ ) expected fol- 
lowing standard cosmological recombination liSeager et alJll99fj) . 
this is a reasonable assumption. Howev er, once aij > 0.01, electron 
collis ions become quite important ('e.g.. lNusseil2005tlKuhlen et alJ 
l2006h . In this section, we will apply our improved calculation of 
Kio to some simple cosmological examples. We will only consider 
the regime Tk ^ 3000 K, where the effects of electronic excita- 
tions of hydrogen are much less than one percent and the Lya 
background from collisional excitations can be ignored (assuming a 
Maxwell-Boltzmann electron distribution and hence neglecting any 
X-ray background). These examples are therefore not particularly 
realistic, but they isolate the major effects of electron collisions and 
so are useful model problems. 

Equation Q shows that perturbations in the density, temper- 
ature, ionized fraction, and velocity all source fluctuations in the 
brightness temperature. Because, to linear order in fc-space, veloc- 
ity perturbations are simply proportional to density perturbations, 
we can write the Fourier transform of the fractional 21 cm bright- 
ness temperature perturbation as ( Furlane tto et al. 2006.') 



521 (k) = {P + + + Pt5j 



(24) 



where 5, 5^, and 5t are the Fourier-space fractional perturbations 
in density, neutral fraction, and Tk, respectively, and ^ is the co- 
sine of the angle between the line of sight and the wavevector k. 
The pi factors are linear expan sion coefficients. For simplicit y, we 
will assume that 5t ~ Q (see lPritchard & Furla netto' 2006a for a 
detailed discussion of temperature fluctuations). The relevant ex- 



Figure 5. Ratio of the 2 1 cm brightness temperature fluctuation amplitude 
to its value in a purely neutral Universe, as a function of Xi. The thin and 
thick curves are for z = 12 and z = 20, respectively. The solid and dashed 
lines take Tk = 100 and 3000 K. The dotted curve shows the fluctuations 
at z = 12 and Tk = 3000 K when 5a: = 0. 



pansion coefficients Pi are then iFurlanetto et alj2006h 

1 



P 



1 + 



(25) 
(26) 



We assume that 5x is determined by photoionization equilib- 



num. 



U +XHI/ 



5 = gxSx 



(27) 



where xm is the global mean neutral fraction. This is not 
a particularly good assumpt ion during reionization (again, see 
IPritchard & Furlanettdl200'63 for a more careful tre atment), but i t 
allows us to compare our work with previous results usse j200^ . 
In any case, 5x generally makes only a small contribution to the 
fluctuations. 

From equation <24t . it then follows that the spherically- 
averaged power spectrum of 21 cm fluctuations can be written 

P2i(fc) = STliP'^ + 2/373 + l/5)Pss{k), (28) 

where 5Tb is the mean brightness temperature, Pss (fe) is the matter 
power spectrum, and P' = P + PxQx- We will quote our results in 
terms of the mean temperature fluctuation as a function of scale, 
Aii(fc) = (fcV2^')P2i(fc). 

Figure |4| shows the brightness temperature fluctuations in 
some example scenarios at z — 12.^ The lower thick solid curve 
shows A21 for the standard calculation, in which the only heat 
source is Compton scattering of CMB photons (computed with 

Note that here we use the linear version of the matter power spectrum, 
PgS, with the transfer function of Eisenstein & Hu ( 1998). At the smallest 
scales in the plot, nonlinear corrections will actually be significant (e.g., 
Illiev et al.2003.:.Furlanetto et al..20Q4:.Naoz & Barkana,2Q05i) . 
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RECFAST: fSeager et alJl999^ . In this case the gas is cold {Tk = 
3.51 K) and almost entirely neutral {Xi = 1.8 x 10~*), so both col- 
lisional coupling and the 21 cm fluctuations are extremely weak. 
The upper thick solid curve shows a contrasting case in which 
Ts = 10^ K and Xi — 0. Here the IGM has saturated in emis- 
sion. 

The thin curves assume Tk — 10^ K and take Xi — 
0, 0.01, 0.05, 0.1, and 0.2, from bottom to top. Naively, of course, 
one would expect A21 to decrease as Xi increases, simply because 
there is less neutral gas. But, as we have seen, H-e~ collisions are 
much more efficient at changing Ts than H-H interactions, so they 
can enhance the fluctuation amplitude by more than a factor of two 
even at these relatively small amounts of ionization. As a result, 
warm gas can seed milliKelvin fluctuations even in the absence of 
a Lya background (regardless of the heat source). 

Figure|5|shows how the amplification factor depends on Xi by 
comparing to the fluctuation amplitude in fully neutral gas. First 
consider the two thick curves, which show the ratio between the 
signals with and without ionization at z = 20. The solid and dashed 
curves assume Tk ~ 100 and 3000 K, respectively. The amplifica- 
tion is larger for cooler gas because the overall coupling is weaker 
at lower temperatures, and the addition of free electrons makes 
more of a difference. Interestingly, electron coupling provides a 
large enough boost that ionizing the gas continues to increase A21 
until Xi is quite large - even as high as 90% at the lower tempera- 
ture. 

The thin curves show the same ratios at 2 = 12. Here the 
amplification is even larger because the densities are smaller (and 
hence the coupling weaker for a given Tk and Xi). For Tk ~ 
100 K, we also show the ratio assuming Sx ~ (rather than deter- 
mined by photoionization equilibrium) with the dotted curve. Ob- 
viously, variations in the ionized fraction play only a minor role 
compared to the density fluctuations. 

Our results can be compared directly with those of iNusseJ 
yOOS), who presented more detailed calculations of 21 cm fluc- 
tuations in a warm Universe (in his case, one flo oded with X - 
ray photons). He used the kio values cal culated bvlFiel j jl958t), 
follow ing the approximate phase shifts of iMassev & MoiseiwitchI 
il95lh . These overestimated the scattering lengths by ^ 25%, so 
our predicted amplitudes are somewhat smaller than his. Note again 
that, because we have neglected fast photoelectrons from X-rays 
jChen & Miralda-Escude|l200d IChuzhov et alJ[2006l). as well as 
temperature fluctuations 'Pritc hard & Furl anetto 2006aj. these Fig- 
ures are probably too simplistic to offer anything more than basic 
intuition. 



5 DISCUSSION 

We have computed the spin de-excitation rates for hydrogen in the 
ground state during electron-hydrogen collisions. Previous calcu- 
lations assumed that spin exchange through elastic S-wave scat- 
tering was the only relevant mechanism. Our analysis showed that 
while spin exchange does dominate spin de-excitation, higher par- 
tial waves contribute significantly to the total de-excitation rate. 
We used newer calculations of t he elastic scat tering phase shifts 
for the L = 0-3 partial waves iWang & Cal laway 1994) to re- 
compute Kio. Our main results are presented in Table Q They re- 
main accurate up to Tk ^ 1.5 x 10* K, where collisions above 
the n = 2 threshold begin to become important at the percent 
level. Our results differ from iSmith (.1966i) and especially from the 



widely-used fit in iLiszd i200lh by several percent over the range 
1 K < Tx < 10* K. 

We have also shown that spin coupling produced directly 
through H-e^ collisions will become a secondary effect at 
Tk ^ 6200 K, because Lya photons produced through collisional 
excitations easily dominate the spin coupling in this regime. Al- 
though the excitation rate to the 2p level (and to other levels that 
cascade through Lya) can be much smaller than the collisional 
spin de-excitation rate, each Lya photon scatters ~ 10^ times be- 
fore redshifting out of resonance. This dramatically boosts the ef- 
ficiency of the radiation background relative to collisions. In prac- 
tice, the Lya background may be important at even lower tempera- 
tures. In our calculations, we have assumed a Maxwell-Boltzmann 
distribution of electron velocities, but even a small population of 
nonthermal high-energy electrons can induce significant Lya cou- 
pling. For example. X-rays produce fast secondary electrons able 
to excite hy drogen collisionally, which can easily produce such a 
background <Chen & Miralda-Escudel2006l:IChuzhov et all200^ . 

We also considered some simple applications of our results 
to the high-redshift IGM. Electron-hydrogen collisions are likely 
to be most important if the radiation background is dominated by 
X-rays, so that the IGM becomes warm and weakly ionized. Such 
scenarios are particularly common when black hole accretion pro- 
vides a significant fraction of the ionizing fl ux iRicotti & OstrikeJ 
'2004»: lRicotti et a l."2005':'Kuhlen et al.'200^. We have seen that the 
typical 21 cm fluctuation amplitude can reach the milliKelvin level 
even without a significant ultraviolet background. Our rate coeffi- 
cients will be useful in assessing the implications of such scenarios 
for the 21 cm sky. 

SRF thanks the Tapir group at Caltech for their hospitality 
while much of this work was completed. This publication has been 
approved for release as LA-UR-06-5395. Los Alamos National 
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operated by Los Alamos National Security, LLC, for the National 
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